Spatio-Temporal Development of Vegetation Carbon Sinks and Sources in the Arid Region of Northwest China

Drylands, which account for 41% of Earth’s land surface and are home to more than two billion people, play an important role in the global carbon balance. This study analyzes the spatio-temporal patterns of vegetation carbon sinks and sources in the arid region of northwest China (NWC), using the net ecosystem production (NEP) through the Carnegie–Ames–Stanford approach (CASA). It quantitatively evaluates regional ecological security over a 20-year period (2000–2020) via a remote sensing ecological index (RSEI) and other ecological indexes, such as the Normalized Difference Vegetation Index (NDVI), fraction of vegetation cover (FVC), net primary productivity (NPP), and land use. The results show that the annual average carbon capacity of vegetation in NWC changed from carbon sources to carbon sinks, and the vegetation NEP increased at a rate of 1.98 gC m−2 yr−1 from 2000 to 2020. Spatially, the annual NEP in northern Xinjiang (NXJ), southern Xinjiang (SXJ) and Hexi Corridor (HX) increased at even faster rates of 2.11, 2.22, and 1.98 gC m−2 yr−1, respectively. Obvious geographically heterogeneous distributions and changes occurred in vegetation carbon sinks and carbon sources. Some 65.78% of the vegetation areas in NWC were carbon sources during 2000–2020, which were concentrated in the plains, and SXJ, the majority carbon sink areas are located in the mountains. The vegetation NEP in the plains exhibited a positive trend (1.21 gC m−2 yr−1) during 2000–2020, but this speed has slowed since 2010. The vegetation NEP in the mountain exhibited only intermittent changes (2.55 gC m−2 yr−1) during 2000–2020; it exhibited a negative trend during 2000–2010, but this trend has reversed strongly since 2010. The entire ecological security of NWC was enhanced during the study period. Specifically, the RSEI increased from 0.34 to 0.49, the NDVI increased by 0.03 (17.65%), the FVC expanded by 19.56%, and the NPP increased by 27.44%. Recent positive trends in NDVI, FVC and NPP have enhanced the capacity of vegetation carbon sinks, and improved the eco-environment of NWC. The scientific outcomes of this study are of great importance for maintaining ecological stability and sustainable economic development along China’s Silk Road Economic Belt.


Introduction
The terrestrial ecosystem is a critical component of the global carbon cycle [1]. It is the main platform for atmospheric CO 2 to enter the terrestrial sphere, and as such is highly vulnerable to climate change and human activities [2,3]. When the vegetation carbon stock is greater than the carbon emitted through soil microbial respiration it behaves as the carbon sink. Whereas the carbon emitted through soil microbial respiration is greater than the vegetation carbon stock and behaves as the carbon source [4]. The net ecosystem productivity (NEP) represents the net carbon exchange between the atmosphere whereas this trend of wind speed within 1-3 m s −1 reversed from 1993 to 2015 [33]. Li et al. [18] detected that the drylands in NWC experienced notable land degradation or vegetation brownness during 1980-2015, and a drying trend has been concluded since the 2000s [25,34], as well as notably enhanced evaporation. Scholars the world over have focused on the region's climate-and human activity-related challenges, searching for viable solutions. Widely regarded as an important area in the global carbon cycle estimation, the sparse vegetation cover, fragile ecosystems of NWC are extreme sensitivity to global climate change. Chen et al. [25] detected that the vegetation coverage in NWC generally exhibited a positive trend during 1982-1998 and then an abrupt decline since 1998. Recent desert-oasis transition zones in the basin of NWC showed a reduction in area in response to the warming climate and human activities [35]. Deng et al. [36] found that the start time for the phenology growing season was 0.66-3.45 d later than that of previous years under drought condition. In addition to damaging ecosystems, the variations in vegetation carbon sinks and sources can exert growing pressures on economic development and social stability. Presently, numerous studies have been conducted on the characteristics of soil organic carbon mineralization and the differences in soil respiration rate under different land use patterns. However, what are the current dynamic changes of the vegetation carbon cycle and ecological security in NWC in the context of climate change and human activities? What are the changes in the ecosystem carbon cycle and the potential causes?
Vegetation carbon cycle in NWC plays an important role in the global carbon balance. Exploring the dynamic changes of carbon sinks and sources in NWC and evaluating the ecological security is crucial to the security of the arid ecosystem, and ecological civilization construction. Therefore, this study has analyzed the dynamic changes of spatial and temporal patterns of carbon sinks/sources in the terrestrial ecosystems of NWC, based on multi-source remote sensing image data, Normalized Difference Vegetation Index (NDVI), land use data and meteorological data from 2000 to 2020. The NEP calculation model and classifications were used to conduct an ecological security assessment and detect the influencing factors. The main objective of this study was to provide a theoretical basis for maintaining both the ecological stability and economic sustainable development of the Silk Road Economic Belt.

Study Area
Northwest China is located in the Eurasian hinterland and covers a total area of 250 × 10 4 km 2 . It sits within the boundaries of 73.55 •  Because of its low levels of precipitation, strong evaporation, and vast distance from the nearest ocean, NWC is one of the driest areas in the world. It has sparse vegetation, which is mainly dominated by mountains and desert. The total oasis area accounts for 9.7%, and the desert and Gobi areas accounts for approximately 80% of the region in NWC. A temperate continental climate dominates the region, bringing annual average temperatures of around 6.58 • C and annual average precipitation of 123.30 mm, and the potential evaporation could be up to 2500-3000 mm. The spatial distribution of NWCs water resources is relatively poor and unevenly allotted. Spatially, the water resources are mainly distributed in the western region, and the water resources are less in the eastern region. The resources themselves are mainly composed of glacier and snow meltwater from the mountains. Drought is a common characteristic of this massive arid zone, which has become more and more vulnerable to climate change. The artificial oasis areas in NWC Int. J. Environ. Res. Public Health 2023, 20, 3608 4 of 23 are less than 10% but carry about 98% of the population and produce 95% of GDP. The oasis areas are the main carrier of human activities and economic and social development in NWC. potential evaporation could be up to 2500-3000 mm. The spatial distribution of NWCs water resources is relatively poor and unevenly allotted. Spatially, the water resources are mainly distributed in the western region, and the water resources are less in the eastern region. The resources themselves are mainly composed of glacier and snow meltwater from the mountains. Drought is a common characteristic of this massive arid zone, which has become more and more vulnerable to climate change. The artificial oasis areas in NWC are less than 10% but carry about 98% of the population and produce 95% of GDP. The oasis areas are the main carrier of human activities and economic and social development in NWC.

Materials
Multiple remote sensing data products from the Moderate Resolution Imaging Spectroradiometer (MODIS) data were used in this study for the period 2000-2020. NDVI data were downloaded from the MOD13Q1, with a spatial resolution of 250 m and a temporal resolution of 16 d. The NPP was calculated using the Carnegie-Ames-Stanford Approach (CASA) model, based on the MCD15A3H, MOD15A2H, MCD12Q1, Terraclimate, GLADAS/T3H datasets. The RSEI adopted in the research was estimated using multiple adopted datasets, including MOD09A1 and MOD11A2. Land use data with spatial resolution of 30 m for the periods 2000, 2010, and 2020 were provided by the Resource and Environment Sciences and Data Center, Chinese Academy of Sciences. The land use dataset includes six categories: cultivated land, industrial land, forest, grassland, water, and unused land. The classification standards were interpreted, classified, and digitized according to the land classification system of the Chinese Academy of Sciences, and the decipherment results were verified in conjunction with fieldwork. In this study, the MODIS products and the climate datasets (TerraClimate, GLDAS) during 2000-2020 were downloaded and processed from the Google Earth Engine (GEE) platform. The NPP and RSEI data were calculated by using the GEE platform, the data were processed and resampled at a spatial resolution of 500 m.
The Shuttle Radar Topography Mission (SRTM V4.1) digital elevation model (DEM) with a spatial resolution of 90 m was provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences. Monthly temperature

Materials
Multiple remote sensing data products from the Moderate Resolution Imaging Spectroradiometer (MODIS) data were used in this study for the period 2000-2020. NDVI data were downloaded from the MOD13Q1, with a spatial resolution of 250 m and a temporal resolution of 16 d. The NPP was calculated using the Carnegie-Ames-Stanford Approach (CASA) model, based on the MCD15A3H, MOD15A2H, MCD12Q1, Terraclimate, GLADAS/T3H datasets. The RSEI adopted in the research was estimated using multiple adopted datasets, including MOD09A1 and MOD11A2. Land use data with spatial resolution of 30 m for the periods 2000, 2010, and 2020 were provided by the Resource and Environment Sciences and Data Center, Chinese Academy of Sciences. The land use dataset includes six categories: cultivated land, industrial land, forest, grassland, water, and unused land. The classification standards were interpreted, classified, and digitized according to the land classification system of the Chinese Academy of Sciences, and the decipherment results were verified in conjunction with fieldwork. In this study, the MODIS products and the climate datasets (TerraClimate, GLDAS) during 2000-2020 were downloaded and processed from the Google Earth Engine (GEE) platform. The NPP and RSEI data were calculated by using the GEE platform, the data were processed and resampled at a spatial resolution of 500 m.
The Shuttle Radar Topography Mission (SRTM V4.1) digital elevation model (DEM) with a spatial resolution of 90 m was provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences. Monthly temperature and precipitation data from the Climatic Research Unit (CRU TSV4.06) dataset product (available in 0.5 • resolution) was utilized to detect changes in precipitation and temperature since 1990.
All of the datasets used in this study are summarized in Table 1. Note: SR (surface reflectance); LST (land surface temperature); FPAR (fraction of photosynthetically active radiation); Pre (precipitation); Tem (temperature); SOL (total solar radiation).

Net Ecosystem Production Algorithm
In this study, the NEP calculation model and classifications were adopted for estimating regional carbon sinks and sources. The scale of vegetation carbon sinks/sources can be expressed as the difference between NPP and soil microbial respiration and carbon emission, without considering other natural environmental and anthropogenic factors. The algorithm to calculate the NEP is: where NEP is the monthly net ecosystem production (gC m −2 ), NPP is the monthly net primary production (gC m −2 ), R H is the monthly soil microbial respiration (gC m −2 ), T is the monthly mean temperature ( • C), and P is the monthly precipitation (mm).

NPP Estimation Using the CASA Model
The vegetation NPP is a key parameter of a terrestrial ecosystem's carbon cycle and energy flow. It can also be used as an important ecological index to quantitatively evaluate the capability of carbon deposition from vegetation. Although some NPP products can be obtained, the limitations for large regions, complex topography, coarse resolution, and continuous observations occur, e.g., from MODIS products which have missing values across the NWC, with some areas entirely unrepresented, particularly in the lower basins.
In this study, the CASA model proposed by Potter et al. [37] was used to calculate the NPP, which was realized on the GEE platform. This model considers the climate condition (i.e., the NDVI and land use data, temperature, precipitation and solar radiation), and the simulation is determined using two variables: absorbed photosynthetically active radiation (APAR) and light energy conversion (ε). Because of its high accuracy, the CASA model has been widely used to estimate regional vegetation NPP. Compared to other products, the simulated NPP from CASA model showed a close match (R = 0.82, p < 0.01) with the NPP from MOD17A3H, which proved better in this study. Detailed information on the model can be obtained from [5,20]. The CASA model is formulated as follows: where NPP (x, t) is the NPP value in pixel x and in month t, APAR is the absorbed photosynthetically active radiation (MJ m −2 ), FPAR is the fraction of photosynthetically active radiation, ε is the light energy conversion (gC M −1 J −1 ), and SOL is the total solar surface radiation (MJ m −2 ).
where FPAR max (=0.95) and FPAR min (=0.001) are independent of the vegetation type, NDVI i,max is the NDVI value corresponding to 95% of the NDVI value, and NDVI i,min is the NDVI value corresponding to 5% of the NDVI value. FPAR and SR are calculated as: where SR i,max and SR i,min are the NDVI i,max and NDVI i,min , respectively.
with α set at 0.5. The light energy conversion ε (gC M −1 J −1 ) was calculated as follows: where T i and T j are the low and high temperature stresses on the efficiency of light use, respectively; w ε is the effects of the water stress; and ε max is the maximum light use efficiency (gC M −1 J −1 ).

Fraction of Vegetation Cover (FVC) Calculation
The FVC is the percentage of the vertical projection of vegetation (including leaves, stems and branches) on the ground area. It is an important indicator of surface vegetation cover status and has a strong positive correlation with the NDVI. In this study, the FVC in NWC during 2000-2020 was calculated using the mixed-pixel dichotomy model. The assumption of this algorithm is that the NDVI value of each pixel consists of vegetation and soil. The FVC is estimated as follows: where FVC is the fraction of vegetation cover (%), NDVI Veg is the NDVI value covered by vegetation pixels, and NDVI Soil is the NDVI value of soil or non-vegetation covered pixels.

Remote Sensing Ecological Index (RSEI) Estimation
RSEI is a method to quantitatively evaluate an ecological environment based on remote sensing information and the integration of multiple ecological factors [38,39]. RSEI indicates that the adjustment of four indicators, namely the greenness index (NDVI), wetness index (WET), heat index (land surface temperature, LST) and dryness (normalized difference bare soil index, NDBSI). RSEI can objectively and quickly evaluate the natural ecological environment by using principal component analysis (PCA) [20], which can effectively eliminate the deviation of weight setting by person. The RSEI value ranges from 0 to 1, with a higher value indicating a better eco-environment condition. In this study, the long-term evaluation of eco-environment changes in NWC during 2000-2020 were evaluated by using RSEI. The RSEI algorithm is calculated as follows: where In addition, the Mann-Kendall trend test [28] was adopted to identify the trends of variables, and the Pearson correlation [29] was used to detect the degree of correlation between the variables.

Spatial Distribution of Vegetation Carbon Sinks/Sources
The spatial distribution of average carbon sinks and sources in NWC is based on NEP classifications. As shown in Figure 2, the majority of the vegetation-covered region in NWC constitutes carbon sources (NEP < 0) and accounts for 65.78% of the total vegetation area ( Figure 2d). Spatially, the NEP values were mainly distributed among mountains, river basins and oases. The high-value areas included the Ili River Valley, the Altai Mountains, the Qilian Mountains, and several river basins, while the low-value areas were mainly distributed in the Taklimakan Desert, Gurbantunggut Desert, Kumtag Desert, and Badain Jaran Desert and its periphery.
The annual average NEP during 2000-2020 ranged from −33.84 to 9.40 gC m −2 yr −1 , giving an annual average of −19.08 gC m −2 yr −1 . The spatial distribution of the annual mean NEP in NWC showed obvious differences ( Figure 2). The annual average NEP in northern Xinjiang (NXJ), southern Xinjiang (SXJ), and the Hexi Corridor (HX) was −2.47, −48.36 and −1.92 g C m −2 yr −1 , respectively. Spatially, 58.46% of the vegetation areas in NXJ were carbon sources, 79.91% of the vegetation areas in SXJ were carbon sources, and 57.47% of the vegetation areas in HX were carbon sources ( Figure 2d). The vegetation NEP in most basins in SXJ exhibited obvious carbon sources, whereas the Irtysh River, Emin River and Ili River in NXJ were carbon sinks, with an annual average NEP of 3.1, 4.9 and 83.1 gC m −2 yr −1 , respectively. −1.92 g C m −2 yr −1 , respectively. Spatially, 58.46% of the vegetation areas in NXJ were carbon sources, 79.91% of the vegetation areas in SXJ were carbon sources, and 57.47% of the vegetation areas in HX were carbon sources ( Figure 2d). The vegetation NEP in most basins in SXJ exhibited obvious carbon sources, whereas the Irtysh River, Emin River and Ili River in NXJ were carbon sinks, with an annual average NEP of 3.1, 4.9 and 83.1 gC m −2 yr −1 , respectively.

Temporal Variations of Vegetation Carbon Sinks/Sources
From 2000 to 2020, the annual average NEP in NWC showed an obvious increasing trend at a rate of 1.98 gC m −2 yr −1 (p < 0.01). As shown in Figure

Spatial Variations of Vegetation Carbon Sinks/Sources in Different Basins
As shown in Figure 2e,f, the majority of the vegetation NEP in NWC during 2000-2020 exhibited a positive trend, especially in the mountains. Some 59.19% of the vegetation NEP in NWC showed an increasing trend (Figure 2g), with the highest increases accounting for 4.33% and 44.29% of the total area. Spatially, more than half of the vegetation NEP (60.61%) in NXJ showed an increasing trend, 5.46% and 41.01% of the area showed a significant or extremely significant increasing trend. In SXJ, 57.17% of the vegetation area experienced a rise in NEP, with significant increases and extremely significant increases accounting for 3.31% and 47.38% of the total vegetation area, respectively. Some 59.82% of the vegetation NEP in HX showed an increasing trend, 3.27% and 47.11% of the area showed a significant or extremely significant increasing trend.

Spatial Variations of Vegetation Carbon Sinks/Sources in Different Basins
As shown in Figure 2e,f, the majority of the vegetation NEP in NWC during 2000-2020 exhibited a positive trend, especially in the mountains. Some 59.19% of the vegetation NEP in NWC showed an increasing trend (Figure 2g), with the highest increases accounting for 4.33% and 44.29% of the total area. Spatially, more than half of the vegetation NEP (60.61%) in NXJ showed an increasing trend, 5.46% and 41.01% of the area showed a significant or extremely significant increasing trend. In SXJ, 57.17% of the vegetation area experienced a rise in NEP, with significant increases and extremely significant increases accounting for 3.31% and 47.38% of the total vegetation area, respectively. Some 59.82% of the vegetation NEP in HX showed an increasing trend, 3.27% and 47.11% of the area showed a significant or extremely significant increasing trend.
From 2000 to 2020, 40.81% of the vegetation NEP in NWC experienced a positive trend. Decreases and extremely significant decreases accounted for 2.02% and 30.14% of the total region area, respectively. Spatially, the decreased vegetation NEP of NXJ accounted for 39.39% of the total area. It was mainly distributed in the junction of the northern Junggar Basin and the Altai Mountains, southern Junggar basin, northern Tienshan Mountains, and the Ili River basin. Some 42.83% of the vegetation NEP in SXJ showed a negative trend, and an extremely decreasing trend was seen in 35.04% of the vegetation NEP. Some 40.18% of the vegetation area in HX showed a negative trend. The vegetation NEP in southeast HX, which includes the Shiyang River basin and the upper reaches of the Heihe River basin, decreased significantly (p < 0.01). Throughout the study period, 28.71% of the vegetation NEP showed an extremely significant decrease.
As shown in Figure 4,

Spatial Distributions and Variations of Carbon Sinks/Sources at Different Altitudes
In this study, the arid region of northwest China is divided into mountainous areas (altitude ≥ 1500 m) and plains areas (altitude < 1500 m). As shown in Figure 4   In terms of the dynamic changes of vegetation NEP in mountain and plain areas du ing 2000-2020 (Figure 6), the NEP in the NXJ mountain region showed the highest increa ing trend, rising at an average rate of 3.17 gC m −2 yr −1 (p < 0.01). In the mountains of SX and HX, the NEP increased at rates of 1.88 and 2.52 gC m −2 yr −1 (p < 0.01), respectively. Th NEP in the plains also increased substantially in NXJ and SXJ (p < 0.01) at a rate of 1.0 In terms of the dynamic changes of vegetation NEP in mountain and plain areas during 2000-2020 (Figure 6), the NEP in the NXJ mountain region showed the highest increasing trend, rising at an average rate of 3.17 gC m −2 yr −1 (p < 0.01). In the mountains of SXJ and HX, the NEP increased at rates of 1.88 and 2.52 gC m −2 yr −1 (p < 0.01), respectively. The NEP in the plains also increased substantially in NXJ and SXJ (p < 0.01) at a rate of 1.06 and 1.71 gC m −2 yr −1 , respectively. In the plain of HX, the increasing rate of vegetation NEP during 2000-2020 was 0.32 gC m −2 yr −1 . The NEP in the plains of NXJ, SXJ and HX all showed slower trends from 2010 onward (especially since 2015). This is in contrast to the NEP in the mountains, which exhibited high increasing rates of 10.92, 5.16 and 5.13 gC m −2 yr −1 , respectively since 2010.

Influence Factors of Vegetation Carbon Sinks/Sources Changes
The vegetation carbon sinks and sources changes in the arid and semi-arid regions are driven by numerous factors, including changes in water, temperature, precipitation, evaporation, radiation, vegetation and human activities. The increased precipitation promoted the vegetation growth in Central Asia in the 1990s, but the warming temperatures Obvious geographically heterogeneous distributions and changes occurred in vegetation carbon sinks and sources. Carbon source areas form the majority and are mainly located in plains, which accounts for 54.49% of the total plain area. During 2020, approximately 54.37% of the vegetation NEP in mountains were carbon sinks, whereas 69.81% of the vegetation NEP in the plains were carbon sources. The proportion of carbon sinks of vegetation in mountains increased from 37.43% to 54.37% during 2000-2020, and the proportion carbon sinks of vegetation in the plains ranged from 19.85% to 30.19% during the same period.

Influence Factors of Vegetation Carbon Sinks/Sources Changes
The vegetation carbon sinks and sources changes in the arid and semi-arid regions are driven by numerous factors, including changes in water, temperature, precipitation, evaporation, radiation, vegetation and human activities. The increased precipitation promoted the vegetation growth in Central Asia in the 1990s, but the warming temperatures and decreased precipitation reversed this trend during 1998-2013 [40]. The "warm-wet" pattern was prevailing in HX with a positive trend in NDVI since 1999 [25]. Affected by the climate changes or human disturbances, the regional vegetation carbon sinks in the arid and semi-arid regions will weaken or strengthen; some regions may also be transformed into carbon sinks or carbon sources, e.g., the conversion of natural vegetation to cultivated land could cause rapid decline in soil carbon concentration. The entire desert-oasis ecotone in the Tarim River basin showed a shrinking trend with a large amount of vegetation translated into arable land or poor land since 1990 [35], which has changed the ecological functions and increased the carbon sources. Under the policy of converting farmland to forest or grassland, the sequestration of soil organic carbon has been enhanced, leading to a fall in the soil greenhouse gas emissions in the Loess Plateau [41].

Climate Changes
Climate changes in arid region interact to impose complex and varying limitations on vegetation activities, leading to more carbon emissions from vegetation, soil, cryosphere and water into the atmosphere [42][43][44]. During a comprehensive interpretation of interactive climatic controls on plant productivity in Central Asia, with drought influencing both the carbon budget and balance [6], Chen et al. [45] and Luo et al. [46] concluded that the precipitation in the mountains of Central Asia was increasing at a greater rate than that of the plains, whereas a drier climate was detected in the desert. The characteristics of vegetation NEP in NWC indicated that the coupling effect of precipitation and temperatures were the main drivers of vegetation growth. Yao et al. [22] found that Xinjiang has experienced significant warming and moistening between 1961 and 2018. Zhang et al. [47] detected that the high mountains in Xinjiang were wetter, whereas the deserts had become drier since the 1980s under a warming climate. Spatially, the climate, vegetation as well as the eco-environment in the desert-oasis transition zone showed more sensitivity to a warming climate; as a result, the Taklimakan Desert has expanded 7.8% [48], and the bare land area decreased at a rate of 0.03% yr −1 since 1999 [25].
In NWC, the soil heterotrophic respiration (RH) in the vegetation has been positively correlated with temperature (R = 0.36, Table 2). During the study period from 2000 to 2020, temperatures in NWC have an increasing rate of 0.14 • C/decade (p < 0.01). Most regions in NWC experienced a warming trend, especially in the far northern portion of NXJ and HX. The rising air temperatures have increased soil respiration and consumed the organic carbon; as a result, it has been consuming the carbon budget of vegetation in NWC. In contrast, the temperatures have shown a rapid rate of increase in Central Asia, which has caused substantial carbon loss [5,6]. Precipitation in NWC showed a rate of 2.28 mm/decade during the study period. The HX region, the south and southeast portions of NXJ, and most areas of SXJ were the main recipients of the increase. Compared to the temperature changes, the vegetation NEP was significantly positively correlated with precipitation, giving correlation coefficients of 0.31. Water is the important limiting factor for ecological stability in arid regions. Hao et al. [49] showed that the vegetation in arid regions is more sensitive to water deficit caused by temperature change due to high water vapor pressure differences and a low soil moisture climate environment. Li et al. [34] found that the vegetation in the Tienshan Mountains of NWC is extremely vulnerable and sensitive to water deficits. Recent studies showed that the temperatures in the mountains are warmer than the lower plains [50,51], which has aggravated the melting of water from glacier and snow [52,53]. Under the influence of climate change, the river runoff in the NWC from the glacier and snow meltwater has greatly increased [52,[54][55][56], and this trend will continue. The increasing river runoff replenishment and regulating effect on river runoff will be further enhanced, providing more water resources for the oasis and natural vegetation, thereby helping to improve the ecological environment.

NDVI and FVC Changes
Vegetation changes have great influences on the vegetation carbon sinks and sources. The recent increased carbon budget in NWC may result from the improved vegetation. Xu et al. [57] detected that the interactive effects of climate change and land use changes help prediction of the persistence of vegetation change. Liang et al. [58] found that 13.56% of oasis land surfaces in NWC during 1975-2010 showed significant positive trends, whereas 6.07% of oasis land surfaces exhibited a significant negative trend. Wang et al. [59] found that most of the vegetation covered areas in HX experienced significant changes during the period 2001−2017. In this study, the NDVI showed an increasing trend (p < 0.01), 88.99% of NDVI exhibited an increase during 2000-2020. The NDVI value rose from 0.17 in 2000 to 0.20 in 2020, representing an increase of 17.65% over 20 years (Figure 7).
that the increased FVC have contributed to the capacity of the carbon budget, as well as increasing the carbon sinks.   Spatially, as shown in Table 3, the higher and low FVC in NWC showed a notable increase a rate of 0.04% yr −1 and 0.05% yr −1 (p < 0.01), whereas a decreasing trend was detected in the high, medium and lower FVC. The high FVC in HX increased by 0.07% yr −1 , whereas the high FVC in NXJ and SXJ decreased by 0.002% yr −1 and 0.03% yr −1 , respectively. Moreover, the low FVC in NXJ, SXJ and HX increased at rates of 0.05% yr −1 , 0.05% yr −1 and 0.06% yr −1 , respectively (p < 0.01). A significant positive correlation was detected between FVC and NEP (R = 0.88, p < 0.01), and between FVC and NPP (R = 0.93, p < 0.01). Both the NDVI and FVC have strong positive trends with NEP, which indicates that the increased FVC have contributed to the capacity of the carbon budget, as well as increasing the carbon sinks.

Human Activities
Human activities affect the carbon sinks and sources in arid and semi-arid regions by changing the types and areas of land use. The impact of human activities on the vegetation NEP in NWC is primarily through the modification of artificial vegetation NEP. Due in large part to climate change and human activities, the vegetation area (cultivated land, forest and grassland) has expanded by 26,531.27 km 2 (3.95%), of which cultivated land has expanded significantly (44.50%) from 2000 to 2020 ( Figure 9). Spatially, the cultivated land areas in NXJ, SXJ and HX have experienced a particularly notable expansion, increasing by 46.01%, 58.55%, and 14.65%, respectively. The fragile eco-environment of NWC is particularly vulnerable to human disturbance. From 2000 to 2020, the forest area decreased significantly (22.80%), shrinking 39.37%, 8.10%, and 0.87% in NXJ, SXJ, and HX, respectively. Meanwhile, the grassland area in NXJ increased by 7.63%, but grassland in SXJ and HX saw a decrease of 3.62% and 0.28%, respectively. Land use changes are another critical factor affecting soil carbon sequestration or decomposition and greenhouse gas emissions [41]. Under a warming climate and the rapid expansion of the irrigated area, the desert-oasis transition zone showed a reduction in area [25]. Aumuti and Luo [60] found that the total oases areas in Xinjiang increased by 35% from 1990 to 2008. Zhang et al. [55] detected that the oases in the Tarim River basin have greatly expanded during 2000−2020. The deforestation, expansion of cultivated land and urbanization have increased the ecosystem feedback to climate on a large scale. It is necessary to design and control the scale of artificial oases development and protect the space of the natural oases.

Ecological Security Assessment
Under global warming and vegetation changes, the ecological security of NWC has also undergone significant changes. Variations of RSEI in NWC during 2000-2020 are shown in Figure 10. As can be seen, the total eco-environment in NWC greatly improved over the 20-year period, and the ecological grade changed from poor to medium. More specifically, the RSEI increased from 0.34 in 2000 to 0.49 in 2020, and its growth rate increased by 0.0016 yr −1 . Spatially, the ecological grades in NXJ ranged from medium to good, with the RSEI increasing from 0.46 to 0.60; the RSEIs overall growth rate reached 0.0012 yr −1 . For SXJ, the eco-environment grade ranged from poor to medium, and the RSEI increased from 0.30 to 0.44, at a growth rate of 0.0022 yr −1 . The eco-environment grades for HX ranged from poor (RSEI = 0.30) to medium (RSEI = 0.44), and the RSEI in the region increased at a rate of 0.0006 yr −1 . By using the CA-Markov model, Hou et al. [61] indicated that the ecological environment in the Tarim River basin is expected to steadily improve from 2019 to 2030.

Ecological Security Assessment
Under global warming and vegetation changes, the ecological security of NWC has also undergone significant changes. Variations of RSEI in NWC during 2000-2020 are shown in Figure 10. As can be seen, the total eco-environment in NWC greatly improved over the 20-year period, and the ecological grade changed from poor to medium. More specifically, the RSEI increased from 0.34 in 2000 to 0.49 in 2020, and its growth rate increased by 0.0016 yr −1 . Spatially, the ecological grades in NXJ ranged from medium to good, with the RSEI increasing from 0.46 to 0.60; the RSEIs overall growth rate reached 0.0012 yr −1 . For SXJ, the eco-environment grade ranged from poor to medium, and the RSEI increased from 0.30 to 0.44, at a growth rate of 0.0022 yr −1 . The eco-environment grades for HX ranged from poor (RSEI = 0.30) to medium (RSEI = 0.44), and the RSEI in the region increased at a rate of 0.0006 yr −1 . By using the CA-Markov model, Hou et al. [61] indicated that the ecological environment in the Tarim River basin is expected to steadily improve from 2019 to 2030.

Ecological Security Assessment
Under global warming and vegetation changes, the ecological security of NWC has also undergone significant changes. Variations of RSEI in NWC during 2000-2020 are shown in Figure 10. As can be seen, the total eco-environment in NWC greatly improved over the 20-year period, and the ecological grade changed from poor to medium. More specifically, the RSEI increased from 0.34 in 2000 to 0.49 in 2020, and its growth rate increased by 0.0016 yr −1 . Spatially, the ecological grades in NXJ ranged from medium to good, with the RSEI increasing from 0.46 to 0.60; the RSEIs overall growth rate reached 0.0012 yr −1 . For SXJ, the eco-environment grade ranged from poor to medium, and the RSEI increased from 0.30 to 0.44, at a growth rate of 0.0022 yr −1 . The eco-environment grades for HX ranged from poor (RSEI = 0.30) to medium (RSEI = 0.44), and the RSEI in the region increased at a rate of 0.0006 yr −1 . By using the CA-Markov model, Hou et al. [61] indicated that the ecological environment in the Tarim River basin is expected to steadily improve from 2019 to 2030.
To further detect the driving factors of NEP in NWC during 2000-2020, the relationship between vegetation index (NDVI, FVC, NPP), temperature, precipitation, RSEI, RH, and NEP were analyzed in this study (Table 2). A close relation was found between vegetation index (NDVI, FVC, NPP) and NEP in the NWC (R > 0.8, p < 0.01), but no obvious correlation was noted in the RSEI, TEM and RH, which indicates that the increase in vegetation NEP of NWC has a close positive-exponential relationship to the improved vegetation.   To further detect the driving factors of NEP in NWC during 2000-2020, the relationship between vegetation index (NDVI, FVC, NPP), temperature, precipitation, RSEI, R H , and NEP were analyzed in this study (Table 2). A close relation was found between vegetation index (NDVI, FVC, NPP) and NEP in the NWC (R > 0.8, p < 0.01), but no obvious correlation was noted in the RSEI, TEM and R H , which indicates that the increase in vegetation NEP of NWC has a close positive-exponential relationship to the improved vegetation.

Conclusions
This study has investigated spatiotemporal changes in the vegetation carbon sinks/sources in the arid region of northwest China and evaluated the ecological security and influence factors of these changes. Our preliminary conclusions are as follows: The distribution of vegetation carbon sinks and sources showed obvious spatial and vertical zonal characteristics. Approximately 65.78% of the vegetation area in NWC were carbon sources. Spatially, 58.46% in NXJ were carbon sources, and 79.91% in SXJ were carbon sources. Of the vegetation area in HX, 57.47% were carbon sources. The carbon sources are mainly distributed in plain areas with an annual average NEP of −47.52 gC m −2 yr −1 , while the carbon sinks are primarily located in mountainous areas with an annual average NEP of 3.53 gC m −2 yr −1 .
The annual average carbon capacity of vegetation in NWC changed from carbon sources to carbon sinks from 2000 to 2020. The annual average vegetation NEP ranged from −33.84 in 2000 to 9.40 gC m −2 yr −1 in 2020, with an increased rate of 1.98 gC m −2 yr −1 . This trend has been more notable since 2010, at a rate of 4.59 gC m −2 yr −1 . Spatially, the annual NEP in NXJ, SXJ and HX rose 2. The capability of vegetation carbon deposition in NWC was enhanced, with the declined carbon sources and enhanced carbon sinks from 2000 to 2020. The spatial distribution of the carbon sinks/sources change matches the improvement in the eco-environment. The entire eco-environment of NWC received a boost, as was evident from the rise in ecological grades from poor to medium. The RSEI increased from 0.34 in 2000 to 0.49 in 2020, resulting in a 72.97% improvement in its eco-environment. There were notable increases in the NDVI (increased from 0.17 to 0.20), the FVC rose from 16.96% to 20.28%, and the NPP jumped by 1.59 gC m −2 yr −1 . Spatially, a better eco-environment has been emerging, especially in the mountainous basins. In the context of climate change and human disturbances, the eco-environment of NWC has an overall positive outlook, with some drawbacks. The spatial distribution of the carbon sinks/sources change matches the spatial distribution of vegetation changes. The NDVI, FVC and NPP also showed improvements, which have contributed to the increase in vegetation carbon sinks and the improvement of the ecoenvironment. Recent improvements in the ecological service function, and the increased ecosystem service value under the climate conditions have a great influence on carbon sinks/sources in NWC.
The variations in vegetation carbon sinks and sources in NWC are controlled by climate change, human activity, topography characteristic, policy, etc. Due to the increasing water resources, natural vegetation and land-type there has been considerable changes. Population increase has stressed the utilization of water resource and promoted the transformation of natural vegetation coverage to artificial and industrial land. There is an urgent need to strengthen the coordinated management of water resources and protection of ecological security in NWC. It is important to protect and restore the ecosystem and maintain the health of the ecological service function. Determine and control the scale of oases to a reasonable spatial area. The recovery of the natural vegetation is needed to enhance the capacity of soil organic carbon and reduce the emissions of soil greenhouse gas. Water resources are the limiting factor for ecological attributes in arid and semi-arid areas in NWC, and they are crucial to support the economic and social development and maintain the stability of the eco-environment. There is a need to strengthen the management and optimize the utilization structure of water resources, and improve the utilization efficiency in the arid region. Optimize water rights management, strengthen oasis land resources and control the development of groundwater. In the process of restoration and improvement of the damaged ecosystem, it is necessary to accelerate ecological restoration and improve ecological functions. Further research is needed to provide a theoretical basis and scientific support for building a life community of mountains, plains, rivers, lakes, forests, cultivated land, grass, sand and ice in harmony between humans and nature.